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Abstract. Cosmic ray antiprotons provide an im- 
portant probe for the study of cosmic-ray propaga- 
tion in tlie interstellar space and to investigate the 
existence of Galactic dark matter. Cosmic rays are 
hampered by the Moon, therefore a deficit of cosmic 
rays in its direction is expected (the so-called Moon 
shadow). The Earth-Moon system acts as a magnetic 
spectrometer. In fact, due to the geomagnetic field the 
center of the Moon shifts westward by an amount 
depending on the primary cosmic ray energy. Paths 
of primary antiprotons are therefore deflected in an 
opposite sense in their way to the Earth. This effect 
aUows, in principle, the search of antiparticles in the 
opposite direction of the observed Moon shadow. 

The ARGO-YBJ experiment, in stable data taking 
since November 2007 with an energy threshold of a 
few hundreds of GeV, is observing the Moon shadow 
with high statistical significance. Using about 1 year 
data, an upper Umit of the p/p flux ratio in the few- 
TeV energy region is set to a few percent with a 
confidence level of 90%. 
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I. Introduction 

The study of cosmic ray (CR) antiprotons may present 
an opportunity to investigate the baryonic asymmetry of 
the universe and to uncover evidence for the existence of 
galactic dark matter. Antiprotons are anyway produced 
by standard nuclear interactions of CR nuclei over the 
interstellar medium (spallation processes). Therefore, the 
observation of p abundance in the CR flux is a key to 
understanding CR propagation and provide information 
complementary to that provided by secondary nuclei 
as Li, Be and B or secondaries of iron. UnUke sec- 
ondary nuclei, antiprotons are tracers primarily of the 
propagation history of the proton, the dominant CR 
component |[l]. Recent measurements of the antiproton 
flux up to about 100 GeV lH, IS, H are consistent 
with the conventional CR model, in which antiprotons 
are secondary particles yielded by the spallation of CR 
nuclei over the interstellar medium. 



Nevertheless, exotic models of primary p production 
not ruled out by low energy measurements yet are avail- 
able. Antiprotons can be produced from primordial black 
hole evaporation or in anti-galaxies. In some scenarios 
models of primary antiprotons production can provide 
a p/p ratio increasing up to the 10% level in the few- 
TeV energy range (see the review paper by Stephen and 
Golden [5|). 

In addition, cosmic ray antiprotons, as well as 
positrons, are considered as prime targets for indirect 
detection of galactic dark matter (see for example |6| and 
reference therein). As an example, some recent analyses 
suggest that the overall PAMELA 111, Q p and e+ 
data and ATIC/PPB-BETS E], S e+ + e" data can 
be reproduced taking into account a heavy DM particle 
(M >10 TeV) that annihilates into W+W' or hh ifTOl . 
This scenario implies that the p/p ratio, consistent with 
the background of secondary production up to about 50 
GeV, increases rapidly reaching the 10^^ level in the 
TeV energy region. 

Linsley [HI and Lloyd-Evans IJ2J in 1985 indepen- 
dently explored the possibility to use the Moon or Sun 
shadows as mass spectrometer in order to measure the 
charge composition of cosmic ray spectrum. In particular 
Linsley first discussed the idea to measure the cosmic 
ray antiprotons abundance exploiting the separation of 
the proton and antiproton shadows. The geomagnetic 
field should deflect the antimatter component of the 
cosmic rays in the opposite direction with respect to 
the matter one. Therefore, if protons are deflected by the 
geomagnetic field towards east, antiprotons are deflected 
towards west. If the energy and the angular resolution are 
respectively low and small enough, we can distinguish, 
in principle, between two shadows, one shifted towards 
west due to the protons and the other shifted towards east 
due to the antiprotons one. At high energy (> 10 TeV) 
the magnetic deflection is too small compared to the an- 
gular resolution and the shadows cannot be disentangled. 
At low energy (wlOO GeV) the well separated shadows 
are washed out by the poor angular resolution, thus 
limiting the sensitivity. Therefore, there is an optimal 
energy window for the measurement of the antiproton 
abundance. 

In this paper we report on the measurement of the 
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p/p ratio in the few-TeV energy region exploiting the 
observation of the galactic cosmic ray Moon shadowing 
effect performed by the ARGO-YBJ experiment. 

II. The ARGO-YBJ experiment 

The ARGO-YBJ detector, located at the YangBaJing 
Cosmic Ray Laboratory (Tibet, RR. China, 4300 m 
a.s.l.), is the only experiment exploiting the full coverage 
approach at very high altitude. The detector is consti- 
tuted by a central carpet ^^74x78 m^, made of a single 
layer of Resistive Plate Chambers (RPCs) with ^92% of 
active area, enclosed by a partially instrumented guard 
ring that extends the detector surface up to ^100x110 
m^. The apparatus has a modular structure, the basic 
data acquisition element being a cluster (5.72x7.64 m^), 
namely a group of 12 RPCs (2.80x1.25 m^ each). 
Each chamber is read by 80 strips of 7x62 cm^ (the 
spatial pixel), logically organized in 10 independent 
pads of 56x62 cm^ representing the time pixel of the 
detector The RPCs are operated in streamer mode with 
a standard gas mixture (Argon 15%, Isobutane 10%, 
TetraFluoroEthane 75%), the High Voltage settled at 7.2 
kV ensures an overall efficiency of about 96% |14|. The 
central carpet contains 130 clusters (hereafter ARGO- 
130) and the full detector is composed of 153 clusters 
for a total active surface of ^--^6700 m^. 

All events giving a number of fired pads Npad > ^trig 
in the central carpet within a time window of 420 ns are 
recorded. The spatial coordinates and the time of any 
fired pad are then used to reconstruct the position of the 
shower core and the arrival direction of the primary. 

The ARGO-YBJ experiment started recording data 
with the whole central carpet in June 2006. The period 
until autumn 2007 has been mainly devoted to instal- 
lation and debugging operations, the duty cycle being 
lower in that period. Since 2007 November the full 
detector is in stable data taking at the multiplicity trigger 
threshold Ntrig >20 and a duty cycle ^ 90%: the trigger 
rate is about 3.6 kHz. 

The reconstruction of shower parameters is split into 
the following steps. First the shower core position is 
derived with the Maximum Likelihood method from the 
lateral density distribution of the secondary particles. 
In the second step, given the shower core position, the 
shower axis is reconstructed by means of an iterative un- 
weighted planar fit being able to reject the time values 
belonging to the non-gaussian tails of the arrival time 
distributions. A conical correction with a slope fixed to 
a = 0.03 rad is applied to the surviving hits in order to 
improve the angular resolution |15|. 

III. Monte Carlo Simulation 

The air showers development in the atmosphere has 
been generated with the CORSIKA v. 6.500 code in- 
cluding the QGSJET-II.03 hadronic interaction model 
for primary energy above 80 GeV and the FLUKA 
code for lower energies lfT6l . Cosmic ray spectra of p. 
He and CNO have been simulated in the energy range 



from 30 GeV to 1 PeV following IHT). The relative 
fractions (in % of the total) after triggering by the 
ARGO-YBJ detector for events with N strip >30 are: 
p~88%, He~10%, CNO'-2%. About 3-10" showers 
have been distributed in the zenith angle interval 0-60 
degrees. The secondary particles have been propagated 
down to a cut-off energy of 1 MeV The experimental 
conditions have been reproduced via a GEANT3-based 
code. The shower core positions have been randomly 
sampled in an energy-dependent area up to lO'^xlO^ 
m^, centered on the detector. Simulated events have been 
generated in the same format used for the experimental 
data and they have been analyzed by using the same 
reconstruction code. 

A detailed Monte Carlo simulation of cosmic rays 
propagation in the Earth-Moon system has been de- 
veloped to compare the observed displacement of the 
shadow with the expectations, so disentangling the effect 
of the geomagnetic field from a possible systematic 
pointing error The algorithm is described in ITSl . 

IV. Data analysis 

For the analysis of the shadowing effect a 10° x 
10° sky map in celestial coordinates (right ascension 
and declination) with 0.1°x0.1° bin size, centered on 
the Moon location, is filled with the detected events. 
The background is evaluated with the so-called time 
swapping method as described in |19|. 

To form conservative estimate of the p/p ratio in 
this preliminary work the analysis refers to about 6-10^ 
events collected only during the stable operation period 
of data taking after the following event selection: (1) 
each event should fire more than 30 strips on the ARGO- 
130 central carpet to avoid any threshold effect; (2) the 
zenith angle of the event arrival direction should be less 
than 50°; (3) the reconstructed core position should be 
inside an area 250x250 around the detector; (4) the 
reduced of the final temporal fit should be less than 
100 ns2. 

A significance map of the Moon shadow is shown in 
Fig. [T| It contains all the events belonging to the lowest 
multiplicity bin investigated (N strip = 30 - 60), collected 
by ARGO-YBJ during the period December 2007 - 
December 2008 (1350 hours on-source). The sensitivity 
of the observation is about 17 standard deviations and the 
Moon is shifted westward by about 1°. This means that 
a potential antiproton signal is expected eastward within 
1.5° from the actual Moon position. The median energy 
of selected events is E50 ~1.4 TeV (mode energy ^ 0.30 
TeV) for proton-induced showers. The corresponding 
angular resolution is ^^1.8° (a detailed analysis of the 
Moon shadowing effect is given in |fT9l ). 

The projection along the North-South (East-West) 
direction of the deficit counts around the Moon is shown 
in the upper (lower) panel of Fig. |2] for Ngtrip = 30 
- 60. The vertical axis reports the events contained in 
an angular slice parallel to the North-South (East- West) 
axis and centered to the observed Moon position. The 
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Fig. 1 : Moon shadow significance map observed by the 
ARGO-YBJ detector in 1350 hours on-source for events 
with N strip = 30 - 60 and zenith angle 6 < 50°. The 
color scale gives the statistical significance. 




Fig. 2: Upper (lower) panel: deficit counts observed 
around the Moon projected along the North-South (East- 
West) axis for N strip = 30-60 (black squares) compared 
to the MC simulation expectations. 



width of this band is ±3.3°. The data are in good 
agreement with the MC simulation and the observed 
shadow is shifted westward by about 1 ° (lower panel), as 
expected. Since at the Yangbajing latitude the effect of 
the geomagnetic field along the North-South direction is 
negligible, a nearly symmetrical and centered projection 
is expected. A detailed analysis of this sort of projections 
and of the energy calibration of the detector is given in 
lfT9l . We stress that the systematic pointing error has 
been taken into account in the upper limit calculations. 

The event and background maps around the position 
of the Moon projected along the East-West axis for 




Fig. 3: The event and background maps around the Moon 
projected along the East- West axis for Nstrip >60 are 
compared. 



any background subtraction. 

V. Results and discussion 

In order to evaluate the p/p ratio, a maximum like- 
lihood fit is performed using the p content as a free 
parameter with following procedure: 

1) the Moon shadow EW projection has been drawn 
both for data and MC. 

2) the MC Moon shadow has been split into a "mat- 
ter" part plus an "antiproton" part, in such a way 
that the total amount of triggered events remains 
unchanged: 

^Mc{fnat) — > ^Mc{r;mat + p) 
= (1 - r)^Mc{mat) + r$A/c(p) 

3) for each matter-to-antiproton ratio, the expected 
Moon shadow EW projection ^Mci^'jiTiat + p) 
is compared with the experimental one via the 
calculation of the likelihood function: 

B 

logC{r) ^Y.^ME^{r)]~E,{r)-ln{m) (1) 

i=l 

where the N^ is the number of experimental events 
included within the i — th bin, while Ei{r) is the 
number of events expected within the same bin. 
The number Ei (r) is calculated by adding the con- 
tribution expected from MC (^Mci^'jiTiat + p)) 
to the measured background. 
This calculation has been performed for two different 
multiplicity bins: Nstrip = 30 - 60 and Nstrip >60. The 
distributions of cosmic ray primary energy contributing 
to these intervals are shown in Fig.|4] From the MC sim- 
ulation it results that the median energy of all selected 
events is E50 = 1.85 TeV for Nstrip = 30 - 60 and E50 
= 4.10 TeV for Nstr^p >60. 

The result of this analysis is shown in Fig. |5] for 



N 



strip 



>60 are compared in Fig. |3] As can be seen, N 



strip = 30 - 60 and one finds r,„i„ = <i>(p)/$(mai) — 
the shadow clearly appears in the event map, without —0.065±0. 078. The result is below a physical boundary 
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Fig. 4: Energy distributions of cosmic ray contributing to 
the multiplicity bins: ^strip = 30 - 60 and Nstrip >60. 




r 



Fig. 5: AlogC as a function of r (the (p) content) for 
Nsfrip = 30 - 60. The dashed line AlogC = 0.5 is used 
to determine the 68.3% central confidence level. 



(the p content must be positive). An upper limit of 
7.4% with 90% confidence level is set using the unified 
approach of Feldmann & Cousin ||201 for Nstrip = 30 
- 60. For higher multiplicities (Ngtrip >60) an upper 
Umit of 7.4% is set with 90% c.l.. Notice that the 
two values are similar, in spite of the difference of 
the multiplicity interval. It is a consequence of the 
combination of the two opposing effects of the angular 
resolution and of the geomagnetic deviation about which 
we said before. With the assumed flux composition 
in the few-TeV energy range of 88% protons and 12% 
heavier nuclei responsible of the observed deficit, these 
limits corresponds to a p/p ratios of 0.09 for Nstrip = 
30 - 60 with 90% c.l. and 0.09 for Nstrzp > 60 with 
90% c.l., respectively. In Fig. |6] the ARGO-YBJ results 
are shown with a compilation of available measurements 

ED. 

VI. Conclusions 

The ARGO-YBJ experiment, in stable data taking 
since November 2007 with an energy threshold of a few 
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Fig. 6: The antiproton to proton flux ratio obtained 
with the ARGO-YBJ experiment compared with all 
available measurements. The solid curve refer to a direct 
production model. The dashed line refers to a model of 
primary p production by antistars |5|, while the dotted 
one refers to the contribution of an heavy DM particle 
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hundreds of GeV, is observing the Moon shadow with 
high statistical significance. Using about 1 year data, a 
preliminary upper limit of the p/p flux ratio is set to 0.09 
with a confidence level of 90% at a median energy of 
1.85 and 4.10 TeV We stress that in the few-TeV range 
this result is among the lowest available. In 3 years of 
data taking the ARGO-YBJ experiment will be able to 
lower this limit down to the percent level, thus excluding 
some of the current antiproton direct production models. 
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